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Abstract 

Attention  is  focussed  on  the  asymptotic  stress  and  deformation  Helds  near  the  edge 
of  a  mode  I  crack  propagating  steadily  in  an  elastic-perfectly  plastic  incompressible 
material  under  plane  strain  conditions.  Features  of  previously  reported  results  on  this 
problem  are  reviewed,  with  a  view  toward  establishing  discriminating  characteristics. 
An  asymptotic  solution  valid  for  all  crack  speeds  is  constructed.  This  solution  has 
the  properties  that  the  angular  variation  of  stress  and  particle  velocity  around  the 
crack  tip  are  continuous  and  the  plastic  strains  are  bounded  at  the  crack  tip.  The 
construction  depends  on  the  boundedness  of  hydrostatic  stress  at  the  crack  tip,  and 
the  conditions  under  which  this  stress  measure  is  bounded  are  discussed  in  an  Appendix. 
The  possibility  of  an  asymptotic  field  with  logarithmically  singular  plastic  strain  at  the 
crack  edge  is  examined.  While  fields  with  this  feature  can  be  constructed  by  admitting 
discontinuities  in  stress  and  particle  velocity,  it  is  shown  that  the  sequence  of  mechanical 
states  experienced  in  the  discontinuities  are  inconsistent  with  the  principle  of  maximum 
plastic  work.  The  dynamic  asymptotic  field  does  not  reduce  to  the  generally  accepted 
quasi-static  asymptotic  field  in  the  limit  as  the  crack  speed  goes  to  zero.  It  has  been 
established  for  the  equivalent  case  of  mode  111  elastic-plastic  crack  growth  that  the 
domain  of  validity  of  the  dynamic  asymptotic  field  vanishes  as  the  crack  speed  vanishes, 
and  the  implication  is  that  the  mode  I  result  has  the  same  property. 


*  Present  address:  Department  of  Mathematics,  Imperial  College,  London. 
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1.  Introduction 


Steady-state  dynamic  propagation  of  a  mode  I  crack  in  an  elastic-perfectly  plastic 
incompressible  material  under  plane  strain  conditions'  is  considered.  The  material  is 
assumed  to  be  governed  by  the  Tresca  yield  condition  and  the  associated  flow  rule.  As¬ 
pects  of  this  problem  have  been  considered  in  the  past  by  Slepyan  (1976),  Achenbach 
and  Dunayevsky  (1981),  Gao  and  Nemat-Nasser  (1983)  and  Lam  and  Freund  (1985). 
Both  Slepyan  (1976)  and  Achenbach  and  Dunayevsky  (1981)  took  elastic  compressibil¬ 
ity  into  account.  They  were  able  to  extract  solutions  valid  very  close  to  the  crack  tip, 
but  their  solutions  were  also  restricted  to  vanishing  crack  speed.  It  should  be  noted 
that  in  both  studies  the  Tresca  yield  condition  was  used  together  with  the  Mises  flow 
rule,  so  that  the  fields  described  are  consistent  with  normality  of  the  plastic  strain  rate 
to  the  yield  surface  only  in  the  limit  of  incompressibility.  In  addition,  both  studies  im¬ 
posed  restrictions  on  the  out-of-plane  deformation  which  arise  from  assuming  normality 
of  the  plastic  strain  rate  to  the  Tresca  yield  surface.  This  results  in  volumetric  plastic 
strains  which  vanish  only  in  the  incompressible  limit.  On  the  other  hand,  Gao  and 
Nemat-Nasser  (1983)  restricted  their  attention  to  the  fully  incompressible  case.  They 
reported  a  solution  with  jumps  in  stress  and  particle  velocity  along  radial  lines  emanat¬ 
ing  from  the  crack  tip  for  all  crack  speeds  between  zero  and  the  elastic  shear  wave  speed 
of  the  material.  Lam  and  Freund  (1985)  obtained  an  extension  of  the  Achenbach  and 
Dunayevsky  (1981)  solution  to  crack  speeds  between  zero  and  the  shear  wave  speed  for 
the  case  of  elastic  incompressibility,  and  examined  the  case  of  elastic  compressibility 
numerically. 

All  previous  investigators  assumed  that  the  limiting  value  of  hydrostatic  stress  at 
the  crack  tip  is  finite.  This  critical  assumption  can  be  used  to  show  that  the  particle 
velocity  in  the  direction  of  crack  propagation,  say  viy  has  the  form  ux  ~  A\w(R/r)  + 
B\[9)  as  r  — »  0  in  terms  of  polar  coordinates  r,  0  centered  at  the  crack  tip.  Slepyan 
(1976)  and  Achenbach  and  Dunayevsky  (1981)  obtained  solutions  which  correspond  to 
the  case  A  =  0,  while  Gao  and  Nemat-Nasser  (1983)  looked  at  the  situation  A  ^  0. 
A  direct  comparison  of  these  solutions  is  possible  only  for  elastic  incompressibility 
and  very  small  crack  speed.  From  a  comparison  on  this  basis  it  is  found  that  the 
solution  of  Gao  and  Nemat-Nasser  (1983)  is  discontinuous,  as  already  noted,  and  that 
the  plastic  strain  components  are  singular  at  the  crack  tip,  while  both  the  solution 
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of  Slepyan  (1976)  and  that  of  Achenbach  and  Dunayevsky  (1981)  include  continuous 
stress  variations  and  bounded  plastic  strains.  For  crack  speeds  up  to  0.1  times  the 
shear  wave  speed,  the  computational  results  of  Lam  and  Freund  (1985)  are  in  good 
agreement  with  the  Achenbach  and  Dunayevsky  (1981)  solution.  Each  of  the  available 
asymptotic  solutions  corresponds  to  plastic  deformation  over  the  full  angular  range 
about  the  crack  tip,  that  is,  there  is  no  elastic  unloading  sector,  and  in  no  case  does 
any  solution  completely  reduce  to  the  corresponding  quasi-static  result  in  the  limit  of 
vanishing  crack  speed.  The  computational  results  of  Lam  and  Freund  (1985),  however, 
suggest  the  presence  of  some  elastic  unloading  downstream  from  the  crack  tip. 

The  problem  is  re-examined  here  in  an  effort  to  resolve  the  discrepancies  among  the 
available  solutions,  or  at  least  to  identify  a  discriminating  characteristic.  It  is  shown  in 
the  Appendix  that,  given  certain  reasonable  assumptions  concerning  the  behavior  of  the 
stresses  and  the  velocities  near  the  crack  tip,  the  hydrostatic  stress  must  be  bounded 
at  the  crack  tip.  As  mentioned  above  (and  shown  in  section  2),  this  implies  that 
the  particle  velocity  in  the  direction  of  crack  propagation  has  a  certain  mathematical 
structure.  The  governing  equations  that  lead  to  the  detailed  features  of  this  structure 
are  obtained  in  section  2.  In  section  3  the  asymptotic  solution  for  the  case  A  =  0  is 
obtained.  This  solution  corresponds  to  the  extension  of  the  Achenbach  and  Dunayevsky 
(1981)  solution  proposed  by  Lam  and  Freund  (1985). 

The  case  A  ^  0  is  examined  in  section  4.  It  is  first  shown  that  a  solution  of  this 
form  satisfying  the  yield  condition  over  the  full  angular  sector  around  the  crack  tip  must 
be  discontinuous ,  as  argued  by  Gao  and  Nemat-Nasser  (1983).  It  is  then  shown  that, 
for  the  specific  case  of  elastic  incompressibility  and  perfect-plasticity,  the  maximum 
plastic  work  inequality  prohibits  the  existence  of  discontinuities  in  the  stress,  strain, 
and  particle  velocity  components.  This  strong  result  indicates  that  the  discontinuous 
solution  of  Gao  and  Nemat-Nasser  (1983)  must  be  ruled  out  if  all  states  experienced 
by  a  material  particle  as  the  crack  tip  moves  by  must  be  consistent  with  the  principle 
of  maximum  plastic  work.  The  elimination  of  discontinuities  leaves  the  introduction 
of  an  elastic  sector  as  the  only  potential  means  of  avoiding  a  region  of  negative  plastic 
work.  It  is  shown  in  section  4,  however,  that  the  introduction  of  such  a  sector  would 
necessarily  violate  the  unloading  condition. 

Because  the  region  of  negative  plastic  work  cannot  be  avoided,  it  must  be  elim- 
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inated.  Inspection  of  the  governing  equations  leads  to  the  conclusion  that  this  may 
be  accomplished  only  by  choosing  A  =  0.  Thus,  the  solution  presented  in  section  3  is 
the  only  available  asymptotic  solution  consistent  with  the  principle  of  maximum  plastic 
work.  In  addition,  the  solution  for  an  elastic  region  obtained  in  section  4  shows  that  for 
A  —  0  any  elastic  region  must  be  a  region  of  constant  (spatially  uniform)  state.  Because 
the  stresses  must  be  continuous  at  an  elastic-plastic  boundary,  any  elastic  region  must 
also  be  at  yield.  However,  regions  of  this  nature  have  been  considered  in  the  solution 
in  section  3,  and  it  is  therefore  concluded  that  this  solution  is  unique  within  the  class 
considered. 


2.  Governing  Equations  in  Plastically  Deforming  Sectors 


Consider  a  mode  I  tensile  crack  growing  steadily  at  speed  u  in  an  elastic-perfectly 
plastic  material  under  plane  strain  conditions.  The  material  is  assumed  to  be  incom¬ 
pressible,  to  yield  according  to  the  Tresca  criterion,  and  to  deform  plastically  according 
to  the  associated  flow  rule.  Both  rectangular  coordinates  £i,£2  and  polar  coordinates 
r,  0  are  introduced  with  their  common  origin  at  the  tip  of  the  crack.  The  coordinate 
systems  translate  with  the  crack  tip  which  moves  steadily  in  the  X\ -direction  with  speed 
u\  see  Fig.  1.  The  equations  of  momentum  balance  are 


drrjj  ^  dv{ 

dxj  P  dt  ’ 


(2.1) 


where  aij  are  the  rectangular  components  of  the  Cauchy  stress  tensor,  Vi  are  the  rect¬ 
angular  components  of  the  particle  velocity,  and  p  is  the  mass  density  of  the  material. 
The  summation  convention  for  repeated  indices  is  adopted,  where  Roman  indices  have 
the  range  1-3.  It  is  convenient  to  introduce  the  dimensionless  quantities 

1  _  1  _  u 

<?ij  =  -Cij  >  vi  =  -Vi  ,  m  ==  —  ,  (2.2) 

p  u  c5 

where  p  is  the  elastic  shear  modulus  of  the  material  and  cs  =  y/p/p .  Under  steady- 
state  conditions  the  material  time  derivative  of  any  field  quantity  /  is  equivalent  to  the 
appropriately  scaled  spatial  derivative 


/  =  —udf/dxi  . 


(2.3) 
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The  momentum  equation  (2.1)  then  becomes 

d(7ij  2  i 

dxj  dxi 


(2.4) 


In  terms  of  the  dimensionless  particle  velocity,  the  rectangular  components  of  the 
small  strain  rate  tensor  are 


_  u  (  dv{  dvj 

2  \dxj  dx{ 


(2.5) 


and  additive  decomposition  of  strain  into  elastic  and  plastic  parts  is  assumed,  that  is, 


*ij  =  , 


(2.6) 


where  and  €?.  are  the  rectangular  components  of  the  elastic  and  plastic  strain  tensors 
respectively.  The  rectangular  components  of  the  elastic  strain  rate,  plastic  strain  rate, 
and  deviatoric  stress  rate  tensors  are  and  ,  respectively.  The  elastic  strain 

components  are  related  to  the  dimensionless  deviatoric  stress  components  ( sij  =  aij  — 
through  Hooke’s  law, 


(2.7) 


The  plastic  strain  rate  components  are  determined  from  the  Prandtl-Reuss  flow 

rule 

t\-  =  Xfxsij  ,  (2.8) 

where  A  is  an  unknown  non-negative  function  of  position.  Combining  (2.5-8)  yields 

u  ( dvi  dvj  \  1  . 

2  +  dbi)  ~ XflSij  +  2Sij  (2,9) 

which  describes  the  material  response  in  plastically  deforming  regions.  The  Tresca 
yield  condition  in  plane  strain  is 

^(^ll  —  ^22  )2  +  <^12  =  k2  (2.10) 

where  k  =  <jy/2jjJ  is  the  normalized  yield  stress  and  <jy  is  the  yield  stress  in  uniaxial 
tension.  This  condition  states  that  in  plastically  deforming  regions  the  maximum  shear 
stress  in  the  plane  of  deformation  has  the  value  cry/2. 
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The  deformation  field  possesses  the  symmetry 


vi(r,0)  =Vl(r,-0)  ,  v2(r,0)  = -v2(r,-0)  (2.11) 


corresponding  to  mode  I  opening.  The  crack  faces  are  traction  free,  which  is  expressed 

by 

022(7’,  ±7r)  =  0  ,  <r12(r,±Tr)  =  0  (2.12) 

for  all  r  >  0. 

As  first  observed  by  Koiter  (1953)  for  elastic  compressibility,  the  assumption  that 
the  out-of-plane  stress  <733  is  the  intermediate  principal  stress  and  the  subsequent  ver¬ 
ification  from  the  solution  implies  that  there  is  no  plastic  strain  in  the  out-of-plane 
direction,  that  is,  e33  =  0.  In  light  of  (2.6-7),  this  implies  that 

=  -^22  ,  (2.13) 


which  allows  the  yield  condition  (2.10)  to  be  written  as 

s22  +  si2  —  k2  •  (2-14) 


This  condition,  with  (2.13),  implies  that  the  deviatoric  stresses  are  bounded  for  all 
values  of  r.  If,  in  addition,  it  is  assumed  that  the  hydrostatic  stress  a  —  l<Tj.fc  is 
bounded,  then  the  momentum  equations  can  be  used  to  discern  the  form  of  tq.  (The 
Appendix  contains  a  discussion  of  the  conditions  under  which  a  may  be  shown  to  be 
bounded.) 


A  typical  momentum  equation  is 

do- 11  d<j\ 2  2  dv\ 

dxi  dx2  m  dxi 


(2.15) 


Tieating  the  dependent  variables  as  functions  of  r  and  0 ,  rectangular  derivatives  may 
be  replaced  with  polar  derivatives  by  a  simple  change  of  variables,  for  example, 


df. 

dx3 


=  cos  0 


dl 

dr 


sin  0  df 
r  dO 


(2.16) 
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If  the  stresses  are  bounded  as  r  — >  0,  then  in  view  of  (2.16),  the  lefthand  side  of  (2.15) 
is  0(l/r).  The  right  side  of  (2.15)  must  therefore  be  0(l/r).  Now  (2.16)  may  be  used 
to  determine  what  form  of  V{  gives  rise  to  the  strongest  allowable  singularity.  If  the 
particle  velocity  is  logarithmically  singular  as  r  — >  0,  then  the  first  term  in  (2.16)  gives 
rise  to  a  1/r  singularity.  In  order  to  avoid  a  stronger  singularity  arising  from  the  second 
term  in  (2.16),  the  coefficient  of  the  lnr  term  must  be  independent  of  9.  Alternatively, 
if  the  particle  velocity  is  a  function  of  9  as  r  — >  0,  then  the  second  term  of  (2.16)  gives 
rise  to  a  1/r  singularity  while  the  first  term  is  o(l/r).  In  summary,  it  appears  that  the 
two  possibilities  leading  to  the  strongest  allowable  singularity  in  particle  acceleration 
are 

.  .  [R\  dAi 

Vi  ~  Ai  In  (- J  ,  —  =  0  (2.17a) 

or 

Vi  -  Bi{9)  (2.176) 


where  R  in  (2,17a)  is  a  parameter  with  the  dimensions  of  length.  R  is  assumed  to  be 
the  same  order  of  magnitude  as  the  maximum  extent  of  the  active  plastic  zone.  The 
plane  strain  condition  implies  that  f?3($)  =  0  and  ^43  =  0. 

It  is  assumed  that  the  behavior  of  may  change  in  different  angular  sectors 
about  the  crack  tip.  Symmetry  requires  that  A2  =  0  in  the  upstream-most  sector,  and 
therefore  the  result  from  section  4.2,  that  the  velocity  components  must  be  continuous, 
implies  that  there  can  be  no  logarithmic  behavior  of  V2  in  any  angular  sector  about  the 
crack  tip.  Thus  (2.176)  is  the  asymptotic  representation  of  V2  for  0  <  6  <  7r,  and  from 
here  on  A\  is  replaced  by  A. 


The  incompressibility  condition  implies  that  the  particle  velocity  is  divergence  free, 
which  in  plane  strain  may  be  expressed  by 


dv1  dv2 
dx1  dx2 


(2.18) 


This  implies  that  the  velocity 
according  to 


components  are  derivable  from  a  velocity  potential  </>(r,  6) 


d<f> 


d(j) 


dx2  ’ 


9xi 


Vl  = 


V2  =  - 


(2.19) 


If  the  velocity  components  are  taken  to  be  a  combination  of  the  forms  indicated  in 
(2-17), 


vi  =  Aln(^)  +Bi{0)  ,  (2.20a) 

V2  =  B2(6)  ,  (2.20 b) 


then  the  divergence  condition  (2.18)  requires  that  A,  i?i($),  and  B2(6)  satisfy  the 
relation 

B[(0)  tan  9  =  B'2(0)  -  A  ,  (2.21) 

where  the  prime  denotes  differentiation  with  respect  to  6.  This  relation  allows  integra¬ 
tion  of  (2.19)  to  obtain 

9)  =  Ax2  (in  (*)  +  l)  +  x2B1(0)  -  x,B2{e)  .  (2.22) 


to  within  an  arbitrary  constant. 

It  should  be  pointed  out  that  the  asymptotic  form  of  Vi  is  actually 

vi  ~  Ain  +  B(r)  +  ^i(^)  +  o(l)  as  r  — >  0  (2.23) 

where  F{r )  =  o(lnr)  does  not  contribute  to  the  leading  order  terms  in  any  of  the  gov¬ 
erning  equations,  and  therefore  cannot  be  determined  without  a  higher  order  analysis. 
(See  for  example  Sirovich  (1971)  for  a  discussion  of  the  Landau  symbols  ’O’,  and  ’o’.) 
The  function  F(r)  has  been  neglected  in  this  paper.  If  such  a  function  were  part  of 
the  correct  asymptotic  form  of  v a  it  would  not  affect  the  results  presented  here  for  the 
stresses,  strain  rates,  or  V2 •  However,  the  results  for  the  strains  and  V\  would  require 
important  modification.  That  there  is  no  similar  function  in  the  asymptotic  form  of 
v2  follows  from  the  argument  above,  which  led  to  the  exclusion  of  a  logarithmic  term 
from  v2  ■ 
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The  material  strain  rates  are  obtained  from  (2.5)  and  (2.20)  as 


en  =  sin 26B'2(0) 

2x  2 

^22  =  — cn  (2.24) 

e12  =  (cos 2^(0)  -  A)  . 

2X2 

The  strain  rate  components  are  therefore  singular  like  1/r  as  r  — >  0. 

In  any  region  of  the  body  that  is  at  yield,  the  yield  condition  (2.14)  can  be  satisfied 
identically  if  the  two  deviatoric  stress  components  appearing  in  (2.14)  are  expressed  in 
terms  of  a  stress  function  0(0)  by 

522  ^  k  cos  (20  —  0(0))  5  S12  =  — fcsin(20  —  V>(0))  •  (2.25) 

The  deviatoric  stress  components  are  defined  in  terms  of  this  particular  function  0(0) 
in  order  to  keep  the  notation  similar  to  that  used  by  Gao  and  Nemat-Nasser  (1983). 
Because  the  hydrostatic  stress  a  is  bounded  (given  the  assumptions  of  the  Appendix), 
it  is  taken  to  be  a  function  of  0  where 

<r(0)  =  +cr22)  (2.26) 

In  terms  of  the  functions  B2(0 ),  -0(0),  and  cr(0),  the  symmetry  conditions  (2.11)  and 
the  boundary  conditions  (2.12)  are 

0(0)  =  0  ,  B2( 0)  -  0  (2.27a) 

and 

0(tt)  =  7 r  ,  <r(7r)  =  k  (2.276) 

respectively,  for  the  range  0  <  6  <  it.  The  system  of  ordinary  differential  equations 
that  these  functions  must  satisfy  is  obtained  next. 

The  governing  equations  that  remain  to  be  satisfied  are  the  two  momentum  balance 
equations  (2.4)  and  the  two  rate  equations  (2.9)  describing  the  material  model.  If  the 
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lepresentations  (2.20),  (2.25),  and  (2.26)  are  substituted  into  these  four  equations  then 
a  lengthy  but  straightforward  manipulation  of  the  results  leads  to 


(2.28) 


B[(9)  tan  9  =  Bl2(0)  —  A 


These  differential  equations  are  considered  for  the  case  (2.176),  that  is  A  =  0,  in  the 
following  section.  It  is  shown  in  section  4  that  it  is  not  possible  to  construct  a  solution 


for  A  7^  0  if  the  principle  of  maximum  plastic  work  is  to  be  enforced. 


3.  Solution  for  A  ==  0 

A  solution  to  the  system  of  equations  (2,28)  is  sought  for  the  case  when  A  —  0. 
This  condition  is  equivalent  to  the  assumption  that  the  particle  velocity  components 
have  the  asymptotic  form  (2.176).  A  solution  of  this  form  has  been  suggested  by  Lam 
and  Freund  (1985).  Attention  is  first  focused  on  the  function  ift(9)  which,  according  to 
(2.28),  must  satisfy 


(i>\9)  -  2)  (cos2  jj)(9)  -  m 2  sin2  6)  —  0 


(3.1) 


m  sin 


in  the  interval  0  <  0  <  tt,  varying  from  V>(0)  =  0  to  ^(tt)  =  tt.  Clearly  (3.1)  is  satisfied 
if  either  0 )  =  2  or  cos  if>(0)  =  ±rasin0.  The  curves  defined  by  cos  i/>(0)  =  ±msind 
are  shown  in  Fig.  2.  Any  line  in  this  plane  with  a  slope  of  2  is  an  integral  curve  of 

W)  =  2. 

A  complete  solution  is  constructed  in  the  following  way.  The  curve  ip(0)  =  20 
satisfies  the  boundary  condition  at  0  =  0,  and  the  curve  =  20  -  if  satisfies  the 

boundary  condition  at  0  —  7r.  The  two  curves  can  be  joined  by  either  of  cosifi(O)  = 
±7?r  sin  0,  and  the  choice  is  governed  by  the  requirement  that  the  plastic  flow  factor 
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\(9)  given  by  (2.28)  must  be  non-negative.  The  factor  (t/>'( 9 )  —  2)  is  negative  along 
both  possible  joining  curves.  Furthermore,  tan  ^(9)  sin  9  =  sin  ifi(d)  where  ~  sin^ 
is  non-negative  along  either  curve.  Therefore,  the  curve  cos  'iji(O)  =  —rosin#  must  be 
chosen  to  render  A  non-negative.  For  the  case  m  =  |,  the  complete  solution  for  ifi(9)  is 
shown  as  the  solid  line  in  Fig.  2,  and  it  is  the  only  continuous  solution  satisfying  the 
boundary  conditions. 


The  interpretation  of  this  solution  is  straight  forward.  The  crack  tip  field  consists 
of  a  uniform  region  in  0  <  0  <  a  non-uniform  region  in  9\  <  9  <  0J*  and  another 
uniform  region  in  9%  <  9  <  7r.  That  the  upstream  and  downstream  regions  are  uniform 
follows  directly  from  ^(9)  —  2  =  0  and  (2.25).  The  transition  angles  are  determined  to 
be 


9\  —  sin" 


m  +  y/S  +  m 2 


$2  =  7T  —  sin  1 


-m  +  V8  +  m2 


(3.2) 


and  the  deviatoric  stress  components  are  determined  by  means  of  (2.25)  from 

(20  ,  i{0<6<6* 

ip(0)  =  <  cos-1  (— m sin 6)  ,  if  6±  <  6  <  6^  (3-3) 

(  26  —  7T  ,  if  $2  <  6  <  TT 


With  ift(6)  determined  the  remaining  differential  equations  in  (2.28)  may  integrated. 
The  first  integrals  are  constant  in  0  <  6  <  6\  and  in  <  9  <  n.  In  the  intermediate 
region,  after  application  of  the  boundary  conditions,  it  is  found  that 


k 

B\{9)  =  (/?2 F(6-,  m)  —  E{6\  m)  +  2m  sin  6 

—fl2F(6\ ;  m)  +  E(y6\ ;  m)  —  2m  sin  )  +  5° 

82(6)  =  (V 1  —  m2  sin  6  —  2m  cos  6  —  y/l~—  m2  sindj  +  2m  cos  9\  j 
a(6)  =  k(  1  +  msind  -  2E(6;m)  -  m sin 6*2  +2E(6*-m)) 


where  F(6]m)  and  E(6-,m)  are  elliptic  integrals  of  the  first  and  second  kind,  respec¬ 
tively,  both  with  parameter  m,  and  (32  =  1  —  m2.  is  a  constant  whose  value  cannot 
be  determined  from  the  asymptotic  analysis. 
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With  the  functions  Bi(0)  determined  by  (3.4),  the  expressions  for  strain  rate  com¬ 
ponents  in  terms  of  B{(9)  (2.24),  may  be  integrated  along  a  line  x 2  =  constant  from 
x1  —  x2  cot  0^  through  decreasing  values  of  x±  to  determine  the  strain  distribution  in 
the  non-uniform  region  9 j  <  9  <  9^.  The  results  of  this  integration,  after  application 
of  the  boundary  conditions,  are 


en(0)  —  —622(0)  —  —Bi(0) 


^12 (0)  = 


m 

"2 


2  (  1  —  \/l  —  m 2  sin  0\  / - r — ; — - 

-  In  1  - : — -  I  +  V 1  —  m2  sm  0  —  2m  cos  0 

\  m  sin  9  I 


m~  ,  /  1  —  ^/l  —  rn2  sin  0^ 


1 

~2~ 


(3.5) 


m  sin  0J 


-  vr  —  m2  sin  0j  +  2 m  cos  0j 


,  ,  0\  /  0* 
-min  tan  -  +  min  tan  — 

2/  V  2 


The  strain  components  are  uniform  in  0  <  0  <  0^,  and  in  9^  <  9  <  ir.  The  uniform 
values  may  be  obtained  from  (3.5). 


4.  Discussion  of  Solution  with  A  ^  0 

I11  this  section,  the  possibility  of  constructing  an  asymptotic  solution  with  loga¬ 
rithmic  dependence  of  Uj  on  the  radial  coordinate  r  is  considered.  If  the  constant  Ay  in 
(2.20a),  is  less  than  zero  the  crack  is  closing.  It  is  therefore  assumed  that  A  is  positive. 

With  reference  to  (2.27-8),  the  stress  function  0(0)  must  satisfy  the  non-linear 
ordinary  differential  equation 

(0;(0)  —  2)  (cos2  ”0(0)  ~  m2  sin2  0)  —  — - Acos(0(0)  —  20)  ,  (4.1) 

rC 


with  boundary  conditions 


■0(0)  =  0 


V’(  7T) 
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=  7T  . 


(4.2) 


The  plastic  flow  factor  A(0)  may  be  manipulated  into  the  form 


2ki 


u 


m 


2  cos  9  cos  ip(9)  +  m2  sin  0  sin (il’(9)  —  20) 
cos2  ip(9)  —  m?  sin2  9 


(4.3) 


Recall  that  the  conditions  (4.1—3)  correspond  to  a  stress  state  satisfying  the  yield 
condition  over  the  full  angular  range  around  the  crack  tip.  This  boundary  value  problem 
for  ip(0)  has  been  considered  by  Gao  and  Nemat-Nasser  (1983),  who  noted  that  a 
continuous  solution  for  ift(O)  is  inadmissible  on  the  grounds  that  it  would  have  to  pass 
through  a  range  of  6  for  which  the  plastic  work  rate  is  negative.  This  fact  may  be 
demonstrated  in  the  following  manner. 

From  (4.3)  it  is  noted  that  the  sign  of  the  flow  parameter  A($)  depends  on  the 
relative  signs  of  the  functions  Aa(0)  and  A2($)  where 


Ai  (#)  =  2  cos  6  cos  ip  +  m2  sin  0  sin(^  —  20)  , 
A2(0)  =  cos2  ip  —  in2  sin2  0  . 


(4.4) 


The  curves  Aj(0)  =  0,  denoted  by  (1,1*),  and  the  curves  A2(0)  =  0,  denoted  by  (U,II' ) 
are  shown  in  Fig.  3;  the  algebraic  sign  of  the  plastic  flow  parameter  is  indicated  in  each 
region. 

Consider  (4.1)  and  the  first  boundary  condition  (4.2).  The  solution  curve  ip(9)  is 
well  defined  up  to  its  point  of  intersection  c  with  the  curve  cos  ip  =  m  sin  9  (curve  II  in 
Fig.  4),  at  which  point  ip'  =  oo.  Equation  (4.1)  can  be  used  to  show  that 

iP'{9)  >  2  ,  ip"(9)  >  0  (4.5) 


for  0  <  9  <  9C,  where  9C  is  the  value  of  9  at  the  point  c.  This  indicates  that  the  solution 
of  (4.1)  up  to  the  point  c  is  of  the  form  shown  by  curve  III  in  Fig.  4,  where  the  line 
ip  =  29  (curve  IV  in  Fig.  4)  is  shown  for  comparison.  In  particular,  (4.5)  shows  that 
the  solution  curve  must  lie  above  the  line  ip  =  29,  at  least  until  point  9  =  9C,  where  it 
encounters  the  curve  cos  ip  =  msh\9  (curve  I  in  Fig.  4)  with  ip' (9)  =  oo.  The  point 
c  will  always  have  the  property  that  ip(9c)  <  7r/2  and,  therefore,  the  segment  of  the 
solution  curve  IE  must  always  lie  inside  the  triangular  region  Odb  defined  by  the  lines 
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0  =  0,  ip  =  7r/2,  and  t/>  =  20.  Point  b  in  Fig.  4  is  the  intersection  of  line  IV  with 
=  7r/2. 

The  only  way  for  the  solution  to  be  extended  continuously  as  a  plastic  state  beyond 
the  point  c  is  through  a  region  of  negative  plastic  work  rate,  which  is  inadmissible. 
Gao  and  Nemat-Nasser  (1983)  proposed  a  remedy  for  this  difficulty  in  the  form  of 
a  discontinuity  in  6 )  initiating  at  the  point  c.  However,  this  approach  is  shown 
later  in  this  section  to  be  inconsistent  with  the  principle  of  maximum  plastic  work,  so 
that  discontinuities  in  field  variables  are  inadmissible  within  the  context  of  the  present 
mechanical  problem. 

Seemingly,  the  only  way  to  construct  a  continuous  asymptotic  solution  with  loga¬ 
rithmic  dependence  of  the  velocity  components  on  r  and  satisfying  the  maximum  plastic 
work  principle  would  involve  the  insertion  of  an  elastically  deforming  sector  starting 
at  some  point  on  the  curve  IE.  This  possibility  is  considered  subsequently  in  section 
4,  where  it  is  shown  that  an  elastic  sector  cannot  possibly  initiate  at  any  point  inside 
the  region  Od6,  and  therefore  on  curve  HI.  The  consequence  is  that  a  solution  with 
logarithmic  behavior  in  the  velocity  components  is  not  admissible,  and  the  velocity 
components  must  therefore  be  bounded. 

Discontinuities 


Possible  discontinuities  in  field  variables  near  the  crack  tip  are  considered  next.  For 
quasistatically  moving  surfaces  of  strong  discontinuity  in  elastic-plastic  solids,  Drugan 
and  Rice  (1984)  have  shown  that  all  stresses  are  continuous,  although  certain  velocity 
components  in  the  plane  of  the  discontinuity  may  suffer  jumps.  The  following  analysis 
shows  that  for  the  dynamic  problem  considered  in  this  paper  all  field  quantities  must 
be  continuous. 

Consider  a  hypothetical  line  discontinuity  £  propagating  with  the  crack  and  in¬ 
clined  at  an  angle  0o  to  the  direction  of  crack  growth.  Define  a  local  rectangular 
coordinate  system  x,y  such  that  y  is  along  the  line  of  discontinuity;  see  Fig.  5.  Let 
the  jump  in  a  quantity  f  be  denoted  by  [/J  =  /+  —  /",  where  the  superscript  plus 
and  minus  indicate  evaluation  of  /  on  the  upstream  and  downstream  sides  of  the  dis¬ 
continuity,  respectively.  Application  of  momentum  balance  across  a  moving  surface  of 
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discontinuity  gives  rise  to  the  following  well  known  jump  conditions 


I <rij\nj  +  m2  —  [v;]  =  0  ,  (4.6) 

where  nj  is  the  normal  to  the  surface  of  discontinuity,  and  c<j  is  the  speed  of  the  surface 
of  discontinuity  in  the  direction  of  the  normal.  From  Fig.  5  Cd  =  usin0o,  and  (4.6) 
may  be  rewritten  as 

KzJ  +  m2  sin0oKl  =  0  ,  Kyi  +  m2  sin0o[uy]  =  0  .  (4.7) 


In  terms  of  the  particle  displacement  components  ux  and  uy,  the  velocity  components 
are 

(  •  /\  O'Uy  dllX  \  /  ,  Oily  Oily  \  .  . 

=  (sm0o  —  -  cos^  -—  J  ,  vy  =  - \*ni0o—  +  cos0o  —  j  ,  (4.8) 

where  use  has  been  made  of  the  incompressibility  condition  dux/dx  +  Ouy/Oy  —  0. 
Continuity  of  displacements  implies,  from  Hadamard’s  lemma,  that 

0  .  (4.9) 


0ux 

Oily 

1  dy  . 

.  &y  1 

Using  (4.9)  in  (4.8)  shows  that 

M  =  0  >  KJ  = 

In  view  of  (4.10),  equations  (4.7)  become 


sin  60 


8Uy 

dx 


-2  sin  #0|eX3/] 


(4.10) 


I0-**!  =  0  >  Kyi  =  2m2  sin2  floKj,] 


(4.11) 


Equations  (2.6-7)  combined  with  (4.11)  imply  that 

K*1  =  0  >  Kyi  =  -4(^,1  ,  Kyi  =2qlepxyj  ,  (4.12) 

where  q  =  ( m 2  sin2  0o)/{  1  —  m2  sin2  0o ). 

It  is  now  assumed  that  material  particles  follow  a  continuous  path  in  strain  space 
through  the  discontinuity  and,  furthermore,  that  if  the  jump  in  a  quantity  is  zero,  then 
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that  quantity  remains  constant  along  the  strain  path  through  the  discontinuity.  In  the 
following  discussion  the  phrase  “ inside  the  jump ”  is  understood  to  mean  ualong  the 

strain  path  from  ej  to  ef.  within  the  discontinuity" .  Thus,  (4.12)  may  be  written  in 
the  differential  form 


daxx  =  0 


dayy  —  ^eyy  > 


daxy  =  2  qdepx 


xy  5 


(4.13) 


inside  the  jump. 


The  maximum  plastic  work  inequality  is 


>0, 


(4.14) 


where  o*j  is  any  stress  state  on  or  inside  the  yield  surface  in  stress  space.  Following 
the  analysis  by  Drugan  and  Rice  (1984)  for  the  quasi-static  case,  the  maximum  plastic 
work  inequality  is  integrated  across  the  discontinuity  to  obtain 


W 


=  JJ  (”•« >0. 


W  is  now  written  in  the  form  W  =  Wp  -  W*  where 

WP  =  Jrl  >  and  w*  =  [l  , 

€ij  J  €* j 

Use  of  (4.13)  allows  evaluation  of  Wp  as 


(4.15) 


(4.16) 


Wp 


=  [**  depxx  +  1  I* 


O' X  y  djO X  y 


1  rjy 

4  U 


^yydcTyy  , 


*»*[«£«]  2q  +  +  g  (aiy  +  Vyy) 


'  2/3/ JJ  J 


(4.17) 


where  the  superscript  S  indicates  the  evaluation  of  a  continuous  field  quantity  at  the 
discontinuity.  Use  of  (4.12)  in  (4.17)  results  in 


w?  = —2(4  +  ■ 


(4.18) 


-16- 


The  integral  W*  is  easily  evaluated  to  give 

w*  =  -<,HI  ■  (4-19) 

Finally,  combining  (4.18)  and  (4.19)  shows  that  W  is  given  by 

W  =  1(2 -  <rj  -  <Ty)l6?f]  >  o  .  (4.20) 

Because  the  material  is  non-hardening,  any  stress  state  in  the  body  may  be  chosen 
for  a*j.  In  particular,  the  choices  a*-  =  crfj,  crT),  or  any  stress  state  inside  the  jump 
are  valid.  Choosing  a*-  =  <r^  and  <j~-  in  (4.20)  gives  >  0  and  [cr^JleUJ  <  0 

respectively.  Clearly  these  relations  are  satisfied  if  and  only  if 

balKl  = 0 

where  use  has  been  made  of  efcfc  —  o.  Using  (4.12)  with  (4.21)  implies  that 

1***1  =  iv?” I*yyl  •  (4.22) 

As  noted  above,  any  state  inside  the  jump  is  a  possible  candidate  for  cr*j  in  equation 
(4.14),  and  therefore  (4.22)  must  hold  for  any  pair  of  states  inside  the  jump.  This 
implies  that 

sxy  =  ±^/qsyy  +  K  (4.23) 

inside  the  jump,  where  K  is  an  unknown  constant.  In  addition,  (4.12)  shows  that  no 
part  of  the  strain  path  through  the  discontinuity  may  be  purely  elastic,  and  hence  the 
yield  condition  must  be  satisfied  at  every  state  inside  the  jump. 

Use  of  (4.23)  in  the  yield  condition  (2.14)  shows  that 

K„J  =  P,„J  =  0  .  (4.24) 

These  results,  together  with  (4.10-13),  show  that  all  field  quantities  must  be  continuous 
across  the  discontinuity,  that  is, 

IK)!  =  0  ,  [Vi]  =  0  .  (4.25) 


(4.21) 
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In  particular,  (4.25)  requires  that  =  0,  and  therefore  that  t/>( 0 )  is  continuous, 

im\  =  o* 


Elastic  Sectors 


In  this  section  the  possibility  of  including  an  elastic  sector  in  the  asymptotic  field, 
initiating  at  some  angle  0 *  on  curve  III  of  Fig.  4,  is  considered.  In  an  elastic  sector 
e\j  =  0  and,  consequently,  the  stress  rate  and  strain  rate  are  related  by 


c ij 


(4.27) 


Use  of  (2.3),  (2.19),  and  (4.27)  in  the  momentum  equations  (2.4)  shows  that  the  velocity 
potential  </>(r,  0)  must  satisfy 

V^)  =  0  (v’=/^  +  fy.  (4.28) 


(Recall  that  f32  —  1  —  m2.)  As  can  be  seen  from  (2.22),  the  velocity  potential  has  the 
general  asymptotic  form 


0(r,  0)  =  rg{0 )  —  Ar  sin  0  In  r  . 


(4.29) 


Representation  (4.29)  implies  that 


V2<£ 


GW 


(4.30) 


where 


G(0)  =  (1  —  to2  sin2  0)  (g(0)  +  gn(0)  —  2A  cos(#)) 


(4.31) 


Using  (4.30)  in  (4.28)  shows  that  G(0)  satisfies  the  ordinary  differential  equation 
G”(0)  +  G(0)  =  0  which,  with  (4.31),  shows  that  g{0)  satisfies  the  inhomogeneous 
equation 

aD  sin  9  +  b0  cos  6  .  , 

g"(0)  +  g(0)  =  - — —  +  2A  cos  0  (4.32) 

1  —  m  sin  d 

where  a0  and  ba  are  constants  to  be  determined.  The  general  solution  of  this  equation 
has  the  form 

g(0)  =  aigi(O)  ,  i  =  1-5  (sum  implied)  (4.33) 
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where  the  Oj  are  undetermined  constants,  and  a3  =  —A.  The  functions  gi(0)  are  given 

by 


(4.34) 


The  strains  may  be  obtained  from  the  velocity  potential  (4.29)  by  means  of  (2.5)  and 
(2.19)  by  integration  along  x2  =  constant.  The  shear  strain  €12  can  only  be  determined 
to  within  an  arbitrary  function  of  x2.  This  unknown  function  is  determined  by  the 
boundary  conditions  at  0  =  9*.  Expressions  for  the  deviatoric  stresses  are  obtained 
through  integration  of  the  constitutive  equations  (4.27), 


2(£ij  -  ««(*»)) 


(4.35) 


and  the  hydrostatic  stress  a  through  integration  of  the  momentum  equations  (2.4). 
Although  the  expression  for  t\2  contains  an  unknown  function  of  x2i  it  turns  out, 
through  use  of  the  momentum  equations,  that  the  difference  e12  -  ep2  in  (4.35)  must 
be  a  constant  multiplied  by  x2.  This  difference  is  asymptotically  negligible,  and  the 
result  is  that  s12  can  be  determined.  The  resulting  expressions  for  the  deviatoric  and 
hydrostatic  stresses  are 


•S11  —  — 2a2  +  2a4  (/?  tan  1  6(j3  tan  0)  —  6)  +  a5  ln(l  —  m2  sin2  6) 
+  2A(l+\nr)-2ep11{x2)  , 


(4.36a) 


(4.366) 


—  m2  a$6  +  m2(A  —  as )  In  r  —  2ep  j  (.t2  ) 


(4.36c) 


where  s0  and  <r0  are  constants. 
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By  means  of  (4.36a.),  boundedness  of  5n  in  the  elastic  sector  requires  that 


eii (^2 )  ~  -A\n(r  sin 0)  +  e0 


(4.37) 


where  e0  is  a  constant.  Also,  because  s12  must  be  bounded,  (4.366)  shows  that  the 
constant  a4  =  0.  Continuity  of  the  hydrostatic  stress  at  an  elastic-plastic  boundary 
implies  that  cr  is  bounded  in  elastic  regions,  and  therefore  as  =  —  A{2  —  m?)/m2 .  The 
resulting  expressions  for  the  stress  components  are 

2  —  m 2 


sn  =  ci  —  A  ^2ln(sin0)  4 - - —  ln(l  —  to2  sin2  0)^  , 


^4 0  — - ~ — tan  J(/?  tan  0)1  , 


s12  ~  C2  2  \  ~  o 

mz  \  p 

a  =  c3  —  2Aln(sin0)  , 


(4.38a) 

(4.386) 

(4.38c) 


where  Ci,  C2,  and  C3  are  constants. 

With  reference  to  Fig.  4,  the  possibility  of  an  elastic  sector  initiating  on  curve  111 
at  some  angle  0  =  0*,  with  0  <  0*  <  0C,  is  considered.  The  purpose  here  is  to  establish 
whether  or  not  it  is  possible  to  avoid  the  region  of  negative  plastic  work  at  0  —  9C 
through  the  introduction  of  a  region  of  elastic  behavior.  In  any  such  elastic  sector  the 
unloading  condition 

—  (^22  +  ^12)  ^  0  (4.39) 

must  hold  in  some  neighborhood  of  0*,  and  in  particular  at  the  point  0  —  0*.  Hence 
the  inequality 

+ *»(**)*«(#•)  ^  °  (4.40) 

must  hold,  where  the  superscript  e  indicates  that  quantities  are  evaluated  on  the  elastic 
side  of  the  elastic-plastic  boundary. 

From  continuity  of  stresses  at  an  elastic-plastic  boundary,  (2.25)  implies 

se22(0*)  =  *  cos  (20*  -  i>*)  ,  se12{0*)  =  -k  sin(20*-^*).  (4.41) 

where  i/>*  =  Use  of  (4.38a),  (4.386),  and  (4.41)  in  the  inequality  (4.40)  shows 
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that  ip*  must  satisfy 


(4.42) 


tan(^*  —  26*)  < 


2  cot  6 *  cos  26* 
4  cos2  6 *  —  m2 


It  can  be  shown  that  this  inequality  is  not  satisfied  anywhere  inside  the  triangle 
Odb  of  Fig.  4.  But,  as  was  mentioned  earlier  in  this  section,  the  point  (ip*,  6*)  must  lie 
on  curve  IE,  which  in  turn  must  lie  inside  triangle  Obd.  Hence  an  elastic  sector  cannot 
be  inserted  into  the  plastic  field  to  circumvent  the  region  of  negative  plastic  work.  This 
implies  that  there  is  no  solution  curve  if  the  region  of  negative  plastic  work  at  6  =  0C, 
indicated  in  Figs.  3-4,  is  present.  Therefore,  ways  to  eliminate  it  must  be  sought. 
Inspection  of  the  governing  ordinary  differential  equations  (2.28)  reveals  that  this  can 
only  be  accomplished  by  choosing  A  =  0,  and  hence  the  velocity  component  rq  may 
not  have  a  logarithmic  singularity  in  r. 

For  A  =  0,  it  can  be  seen  from  (4.38)  that  the  only  possible  elastic  sector  is  a  region 
of  constant  stress  which  must  be  at  yield  in  light  of  full  continuity  at  an  elastic-plastic 
boundary.  But  such  constant  state  regions  are  included  in  the  solution  presented  in 
section  3,  and  it  is  therefore  concluded  that  the  solution  in  section  3  is  unique  within 
the  solution  class  considered. 


5.  Discussion 

There  are  several  reasons  for  the  study  of  the  asymptotic  crack  tip  field  for  dynamic 
growth  of  a  mode  I  crack  in  an  elastic-plastic  material.  To  begin  with,  numerical 
methods,  which  provide  the  only  means  for  obtaining  full  field  solutions  within  this 
problem  class,  are  of  limited  reliability  very  close  to  the  crack  edge.  The  ability  to 
match  computed  fields  to  asymptotic  fields  in  some  sense  establishes  confidence  in  the 
numerical  results.  This  point  was  discussed  by  Lam  and  Freund  (1985).  Secondly, 
the  influence  of  material  inertia  on  the  distribution  of  stress  and  deformation  near  the 
crack  edge  is  of  interest  in  assessing  mechanisms  of  crack  advance  because  these  fields 
represent  the  environment  in  which  the  mechanisms  are  operative.  Thirdly,  in  the 
situation  considered  here,  more  than  one  solution  has  been  suggested  in  the  literature. 
Thus,  an  additional  objective  here  has  been  to  identify  a  discriminating  feature  among 
the  proposed  asymptotic  solutions. 
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The  angular  variations  of  the  stress  components  for  nondimensional  crack  speeds 
m  =  0.1,  and  0.9  are  shown  in  Fig.  6.  The  influence  of  inertia,  aside  from  increasing 
the  transition  angles  0^  and  0 £  somewhat,  is  merely  to  cause  a  moderate  reduction  of 
the  hydrostatic  stress  a  in  the  upstream  region.  Of  course,  this  implies  that  both  <j\\ 
and  (j22  are  similarly  reduced.  Figure  7  shows  the  variation  of  the  strain  components 
for  m  —  0.1,  0.5,  and  0.9.  Inertial  effects  are  more  significant  here  than  for  stress;  the 
downstream  value  of  ej2  is  reduced  almost  to  the  upstream  level,  and  the  variation  of 
both  6x2  and  622  is  reduced  in  the  transition  region.  This  results  in  very  little  strain 
variation  for  large  m.  Note  that  the  unknown  constant  has  been  chosen  to  be 
B°  =  1  for  displaying  e22*  can  on^y  be  determined  by  matching  this  solution  to  the 
outer  solution.  However,  it  might  expected  that  B°  — >  0  as  m  — >  1. 

The  solution  presented  here  includes  no  elastic  unloading  in  any  angular  sector 
about  the  crack  tip,  and  therefore  it  cannot  reduce  to  the  corresponding  quasistatic 
solution  of  Rice,  Drugan,  and  Sham  (1980).  A  somewhat  similar  situation  arises  in 
mode  IH,  where  the  dynamic  steady-state  solution  fails  to  reduce  to  the  corresponding 
quasistatic  solution.  However,  this  issue  was  resolved  by  Douglas  and  Freund  (1982) 
when  they  demonstrated  that  the  domain  of  validity  of  the  dynamic  asymptotic  solution 
vanishes  with  vanishing  crack  tip  speed.  Carrying  this  observation  over  to  the  mode  I 
problem,  it  appears  that  here  too  the  domain  of  validity  of  the  dynamic  asymptotic  field 
vanishes  with  vanishing  crack  tip  speed.  As  pointed  out  by  Lam  and  Freund  (1985), 
the  solution  presented  here  is  an  extension  of  the  Achenbach  and  Dunayevsky  (1981) 
solution  to  the  full  crack  speed  range  0  <  m  <  1  for  elastic  incompressibility.  The 
stress  components  reduce  to  those  observed  in  the  Prandtl  punch  field  as  m  0,  as 
noted  by  Achenbach  and  Dunayevsky  (1981).  The  components  of  strain  (elastic  strain 
is  bounded)  become  unbounded  as  m  0,  although  ei2  remains  zero  in  the  upstream 
sector.  The  region  where  e22  is  singular  as  m  — ►  0  depends  011  the  behavior  of  B%  as 
m  — >  0. 

Concerning  the  previously  proposed  asymptotic  solutions,  it  has  been  noted  that 
the  main  difference  arises  from  whether  or  not  discontinuities  in  stress  and/or  particle 
velocity  components  are  admitted  across  lines  emanating  from  the  crack  tip.  As  is 
shown  in  the  foregoing  discussion,  if  it  is  required  that  the  strain  histories  must  abide  by 
the  principle  of  maximum  plastic  work  through  the  jump,  then  such  discontinuities  must 
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be  ruled  out.  Consequently,  only  asymptotic  solutions  that  exhibit  continuous  angular 
variations  in  stress  and  particle  velocity  are  in  full  compliance  with  the  maximum 
plastic  work  principle. 

The  question  of  whether  or  not  the  strain  paths  through  the  jumps  should  be 
required  to  satisfy  the  principle  of  maximum  plastic  work  could  be  raised,  of  course. 
Within  the  context  of  gas  dynamics,  Courant  and  Friedrichs  (1948)  show  that,  for  a 
mechanical  analysis  of  the  propagation  of  discontinuities,  the  sequence  of  mechanical 
states  experienced  by  a  material  particle  as  a  discontinuity  propagates  across  it  must 
be  the  same  as  if  the  transition  from  the  initial  state  to  the  final  state  has  occurred  in 
a  smooth  simple  wave.  The  reason  for  this  outcome,  in  a  somewhat  simplified  form,  is 
that  in  a  mechanical  description  of  the  phenomenon  (as  opposed  to  a  more  complete 
thermodynamic  description)  the  behavior  of  the  material  does  not  have  the  flexibility 
afforded  by  thermodynamic  properties  to  depart  from  the  mechanical  constitutive  law 
imposed.  These  observations  may  be  carried  over  to  the  present  case  of  elastic-plastic 
crack  growth.  The  model  is  strictly  mechanical  so  it  may  be  expected  that  the  behavior 
of  the  material  in  a  jump  discontinuity  must  also  be  representable  as  the  limit  of  a 
smooth  wave  as  the  slope  of  the  smooth  wave  becomes  indefinitely  steep.  If  this  is  so, 
then  it  follows  immediately  that  the  sequence  of  states  in  the  jump  must  be  consistent 
with  the  principle  of  maximum  plastic  work. 

It  should  be  pointed  out  that  Courant  and  Friedrichs  (1948)  only  state  that  for  a 
mechanical  description  of  the  transition  through  a  discontinuity,  the  initial  and  final 
states  are  the  same  as  if  the  transition  had  occurred  in  a  simple  wave.  However,  their 
proof  of  this  statement  also  demonstrates  that  the  sequence  of  states  in  the  discontinuity 
must  also  be  the  same  as  if  the  transition  had  occurred  in  a  simple  wave.  In  order  to 
describe  the  state  of  a  material  particle  passing  through  the  discontinuity,  they  used 
essentially  the  same  assumption  that  was  used  by  Drugan  and  Rice  (1984),  and  that 
was  used  in  section  4  of  this  paper,  namely,  that  if  the  jump  in  some  quantity  is  zero, 
then  that  quantity  remains  constant  along  the  path  through  the  discontinuity. 
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Appendix:  Bounded  Hydrostatic  Stress 


In  this  section  it  is  shown  that  under  certain  assumptions  on  the  behavior  of  the 
field  variables  and  their  first  two  derivatives  with  respect  to  r,  the  hydrostatic  stress 
a  is  bounded  at  the  crack  tip.  To  this  end,  the  momentum  equations  (2.4)  may  be 
referred  to  crack  tip  polar  coordinates  as 


da 

m2  dv2 

dsrr 

1  d  s-pQ 

dr 

r  86 

dr 

r  86 

-  Z -  5 

r 

1  da 

2  dv2 

1  ds  69 

dsr0 

1 

1 

771  dr 

r  86 

dr 

~  z -  ? 

r 

(A.l) 


where  (2.3)  and  the  incompressibility  condition  have  been  used  to  write 


u 


dv2 


Wr  — 


86 


<A>0 


—u 


:  dv2 

dr 


(A.2) 


and  LOr  and  u)$  are  the  radial  and  circumferential  components,  respectively,  of  tlie 
particle  accelerations.  Differentiation  of  the  second  of  equations  (.A.l)  with  respect  to 
r  gives 


m 
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d2v2 

dr2 


I  (?1 

r2  \de 


ds$Q 

~df 


1  /  d2c^  _  d^see_  dsre  \  d^sre_ 
r  \drdd  drde  +  dr  )  +  dr2  ’  1  j 


while  differentiation  and  manipulation  of  the  momentum  equations  in  rectangular  co¬ 
ordinates  yields 


m2  V2v2 


d2s12  +  d2s  12 


d2 


s  22 


(^.4) 


dx\  dx\  8x18x2 
It  is  assumed  that  all  functions  and  their  first  and  second  derivatives  have  limits  as  r  — » 
0,  and  further  that  all  functions  g(r,  6)  have  the  property  that  dg(r,  6)/ 86  =  0  (g(r,  #)) 
(See  for  example  Sirovich  (1971)  for  a  discussion  of  the  Landau  symbols  ‘O’,  and  co’.) 
Then  (A.l),  (A. 3)  and  (A. 4)  may  be  replaced  by 
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respectively.  A  solution  is  now  sought  of  the  form 


&(r,6)  =  f(r)g(9)  +  h(r,9 )  ,  as  r  — »  0 


(A.8) 


where  f(r)  is  singular  as  r  — *  0  and  /i(r,  0)  =  o(/(r))  for  any  Substituting  (A. 8) 
into  the  second  of  equations  (A. 5)  and  integrating  gives 


m2  v  2 


(r,0)  =  -g\0)  jMdr-J\fedr  +  0^r)  .  (A. 9) 


Sirovich  (1971)  presents  a  discussion  of  the  conditions  under  which  the  order  relations 
may  be  integrated.  Substitution  of  (A. 9)  into  the  first  of  equations  (A. 5)  and  integrating 
once  again  gives 

*M)  =  -9'\9)  J  (~J  ^rdr)  dr~  $  (;  /  lwdr) dr  +  0  (ln2r)  '  (A10) 


Comparison  of  (A. 8)  with  (A. 10)  reveals  that 

g(Q)f{r)  =  -  g"{9)  J  (1  J  ^^dr)  dr 

-J  (;  Jlj^dr^dr  +  0{ln2r)+°(f(r))  ■ 


(A.  11) 


If  g"(9)  7^  0,  then  from  (A. 11)  one  of  the  following  three  conditions  must  hold; 


(*) 


JGJ 


1  f  f(r ) 


dr  )  dr  =  0 


/(:/ 


i  r  i  d2h 
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dr  )  dr 


(«) 


J  {l  J  ~~ dr  =  O  (In2  r)  , 


(*’*’*)  9(d)f(r)  =  -g"(d)  J  (4  J  ~rdr^j  dr  +  °(/(r))  • 

By  previous  assumptions  d2h(r,9)/d92  =  0 (h(r,9))  =  o(f(r))  and  therefore  condition 
(f)  implies  that  J  (l/r  J  f(r)/rdr )  dr  —  o  (f  (l/r  f  f{r)jrdr )  dr)  which  is  clearly  not 
possible.  Furthermore,  successive  division  by  r  and  integration  with  respect  to  r  of 
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the  relation  1  =  o(/(r))  shows  that  In2  r  =  o(/(r)/r).  Thus  (it)  also  implies  that 
J  (l/r  f  f(r)/rdr )  dr  =  o  (f  (l/r  f  f(r)/rdr )  dr)  which  is  impossible.  Condition  (Hi) 
implies  that  f(r )  ^  a  f (l/r  f  f(r)/rdr )  dr,  and  #"(0)  +  0:5(0)  =  Because  neither 
condition  ( i )  nor  (ii)  is  possible  it  is  concluded  that  either 

g"(0)  =  0  ,  or  g"(0)  +  ag(0)  =  0.  (A. 12) 


The  solution  for  <7  is  now  examined  in  the  neighborhood  of  0  =  7 r.  The  constitutive 
equation  (2.9)  shows  that 


dv\  dv2  _  2p.  ^  ds\2 
dx2  ^  dxi  u  12  dxi 


(A.13) 


Differentiation  of  (A. 13)  with  respect  to  X\ 
gives 


d2v2  d2v2  _  d 

dx\  dx2  dxi 

Combining  (A. 14)  with  (A. 7)  shows  that 


and  use  of  the  incompressibility  condition 
^2pAs12  —  •  (^-14) 
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dxi 


u  ds12  \ 
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(-4-15) 


The  traction  free  crack  face  condition  requires  ^(r,  7r)  =  0. 
tion  with  (A. 15),  that 


This  implies,  in  conjunc- 


(4.16) 


Substitution  of  (A. 8)  into  (A. 6)  shows  that 


d_  fldh\ 

dr  \r  dd  ) 


If  it  is  assumed  that 


(4.17) 


(A.1S) 


then  integration  and  multiplication  by  r  implies  that  /(r)  =  o(/(r))  which  is  clearly  not 
possible.  Similarly,  the  assumption  that  (/(r)/r)'  =  0(l/r2)  leads  to  the  conclusion 
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that  f(r)  =  0(1),  which  contradicts  the  assumption  that  f(r)  is  singular.  Therefore 
[A,  17)  may  be  rewritten  as 
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:d2*>2 

8r 2 


=  -9(«)  ( 


f(r)\ 


+  O 


Comparing  (^4.16)  to  (>1.19)  shows  that 


»'WIt),  =  0(?i+o 


(-4-19) 


{A.  20) 


Integration  of  (A. 20)  with  respect  to  r  and  multiplication  by  r  shows  that  either  f(r)  ~ 
0(1)  or  g\i r)  =  0.  Because  f(r)  is  singular  it  must  be  that  g'(7r)  —  0.  In  addition, 
the  traction  free  crack  face  condition  implies  that  <r(r,  7r)  —  —s 22(^5^)-  The  deviatoric 
stresses  are  bounded  for  all  r,  and  therefore  (>1.8)  implies  that  g(7r)  =  0.  Thus,  if  f(r) 
is  singular  then  g(9)  must  satisfy  (A. 12)  with  the  boundary  conditions 

g( 7t)  —  0  ,  and  g{^)  —  0  •  (A. 21) 


The  only  possible  solution  to  either  of  equations  (A. 12),  which  satisfies  (>1.21),  is 
g(9)  =  0.  Thus,  if  /(r)  is  singular  then  g(9)  =  0,  which  implies  that  <r(r,  0)  may  not 
be  singular  in  the  region  near  0  =  it  if  it  can  be  represented  by  (A. 8).  This  result 
is  independent  of  whether  the  material  is  deforming  elastically  or  plastically  in  that 
region.  I11  addition,  it  was  shown  in  section  4  that  the  hydrostatic  stress  must  be 
continuous  throughout  the  body.  Therefore  a  must  be  bounded  in  all  angular  sectors 
about  the  crack  tip. 
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Figure  Captions 


1.  Diagram  of  crack  tip  with  rectangular  and  polar  coordinate  systems.  Both 
coordinate  systems  move  with  the  crack  tip  in  the  x\  direction  with  constant  speed  u. 

2.  Solution  curves  for  ij>(9)  of  the  type  cos  ip(0)  =  d=  sin  9  for  m  =  0.1,  0.5,  and 
0.9.  The  solid  line  represents  the  complete  solution  for  m  =  0.5  and  A  —  0. 

3.  Diagram  of  the  regions  in  which  the  plastic  flow  factor  A  is  positive  (  +  ),  and 
negative  (— ),  for  m  =  0.5  and  d  /  0. 

4.  Curves  which  delineate  regions  where  the  plastic  flow  factor  A  is  either  neg¬ 

ative  or  positive  (see  Figure  3),  and  a  possible  solution  curve  IH.  The  curve  HI  was 
obtained  by  numerically  integrating  the  governing  equation  for  for  m  ~  0.5  and 

.A  =  1.  Curve  IV  is  the  line  ip  =  29. 

5.  Discontinuity  oriented  at  0^  with  respect  to  the  crack  tip,  with  local  x,y 
rectangular  coordinate  system.  The  upstream  (  +  )  and  downstream  (  — )  regions  are  as 
indicated. 

6.  Variation  of  certain  stress  measures  for  0  <  9  <  tt.  All  stress  quantities  have 
been  normalized  by  the  elastic  shear  modulus  fi. 

7.  Variation  of  strain  components  for  0  <  0  <  ix.  Note  that  the  upstream  value 
of  6j i  has  been  taken  to  be  1  for  all  crack  speeds. 
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